Ground state cooling of a nanomechanical resonator in the non-resolved regime via 

quantum interference 
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Ground state cooling of a nanomechanical resonator coupled to a superconducting flux qubit is 
discussed. By inducing quantum interference to cancel unwanted heating excitations, ground state 
cooling becomes possible in the non-resolved regime. The qubit is modelled as a three-level system 
in A configuration, and the driving fluxes are applied such that the qubit absorption spectrum 
exhibits electromagnetically induced transparency, thereby cancelling the unwanted excitations. As 
our scheme allows to apply strong cooling fields, fast and efficient cooling can be achieved. 

PACS numbers: 85.85.-l-j, 07.10.Cm, 85.25.-j, 42.50.Gy 



Nanomechanical resonators (NAMR) currently attract 
considerable interest because of their combination of high 
natural frequencies and large quality factors together 
with a wide range of potential applications Among 
them are measurements of displacement at the quantum 
limit 0], mass measurements biological sensing Q 
and the observation of quantum mechanical phenomena 
in mesoscopic objects To fully utilize the proper- 
ties of NAMRs or to observe mesoscopic quantum phe- 
nomena, it is typically necessary to cool the NAMR to 
the mechanical ground state. Thus it is not surpris- 
ing that a number of different approaches for cooling 
micro- and nanomechanical resonators have been pro- 
posed theoretically and also demonstrated ex- 
perimentally |lll4l7| . For micromechanical resonators, 
cavity-assisted radiation pressure cooling has been in- 
tensely studied 0, 0, [Til, [ll ■ A different approach 
is active feedback cooling, which however typically re- 
quires difficult and precise measurements in real time of 
the displacement of the resonator 

[13, El-El- Cavity- 
based schemes are limited by diffraction, if the size of 
the resonator is small compared to the wavelength of the 
light. For NAMR, it has been proposed to achieve cool- 
ing by periodic coupling to a superconducting qubit (SQ) 
such as a Cooper pair box (CPB) Q or to a three-level 
flux qubit 0. Both techniques rely on a strong reso- 
nant interaction between resonators and the qubit. Re- 
cently, sideband cooling of micro- and nanomechanical 
resonators has attracted considerable interest. For exam- 
ple, cooling a NAMR has been proposed by embedding 
a quantum dot in the resonator [7[, and it was observed 
in a microresonator 11| and in a transmission line res- 
onator [3, El- Also, a quantum theory of cooling has 
been developed [1]. 

A number of problems associated with cooling NAMR 
are shared by laser cooling of atoms or ions. In par- 
ticular, ground state sideband cooling is possible only 
in the resolved regime |8|, in which the motional side- 
bands are resolved from the linewidth of the involved 
transitions Isl . This has been realized recently in 

few systems [la, ll8|, but still this regime typically is 



difficult to achieve, and limits the accessible parameter 
range. To overcome this limit in atomic systems, a cool- 
ing scheme based on electromagnetically induced trans- 
parency (EIT) has been proposed (2l[ and ex- 
perimentally verified in ions 23[. EIT cooling works in 
the non-resolved regime, but suppresses the carrier exci- 
tation without change in the motional quantum number. 
This is achieved by designing the optical properties of the 
target system in such a way that absorption vanishes at 
the carrier transition frequency. 

In this Letter, we discuss ground state cooling of a 
NAMR in the non-resolved regime. The NAMR is em- 
bedded in the loop of a fiux qubit. The qubit is modelled 
as a three-level quantum system in A configuration, and 
time-dependent magnetic fiuxes (TDMF) are applied to 
the qubit in such a way that detrimental carrier excita- 
tions without change in the motional quantum number 
(|n) — >■ \n)) are suppressed by quantum interference. We 
find that the cooling limit of the NAMR has two contri- 
butions. One originates from the scattering of the cooling 
fields, whereas the other one arises from the equilibrium 
phonon number of the environmental thermal bath. Our 
interference-based cooling scheme extends to strong cool- 
ing fields, and thus enables rapid cooling to a high occu- 
pancy of the mechanical ground state. Unlike backaction 
cooling [sl, no significant coherent shift occurs in the fi- 
nal occupancy of resonator in EIT cooling. Our system 
allows for large Lamb-Dicke (LD) parameters via con- 
trolling the applied magnetic field or the working point, 
which leads to rapid cooling to a high occupancy of the 
motional ground state of the resonator. An experimental 
implementation is facilitated by a rather small required 
input power. 

We start with a discussion of the main idea. As shown 
in Fig. [TJa), the NAMR is part of the qubit loop. As it 
vibrates, the area of the fiux qubit changes, which leads 
to a coupling of NAMR and qubit. The qubit is mod- 
elled as a three- level system as shown in Fig. [TJb). The 
driving field fig acts as cooling field, while the field fte is 
a control field. The relevant qubit-resonator level scheme 
in the LD limit is shown in Fig. [ijc). So-called carrier 
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FIG. 1: (Color online) (a) A nanomechanical resonator (red 
curve) coupling to a superconducting flux qubit. (b) The level 
diagram of the flux qubit. (c) Effective level scheme in the 
Lamb-Dicke (LD) limit, (d) Example absorption spectrum of 
the cooling field. Indicated are the transition frequencies for 
the different absorption channels in the LD limit. 



transitions \g,n) — )■ \a,n) do not change the resonator 
state while sideband transitions \g,n) — >■ \a,n± 1) 
do. Both carrier transitions and sideband transitions 
\g,n) — >■ |a,n + 1) with subsequent spontaneous decay 
on average induce heating. In resolved sideband cooling, 
in which the resonator phonon energy hv exceeds the 
transition width h'j, a low cooling limit is achieved since 
the cooling sideband can almost selectively be driven. In 
the non- resolved regime {v < 7), scattering on the car- 
rier transition dominates the heating processes, prohibit- 
ing ground state cooling. The relative probabilities of 
the excitation channels are determined by the absorption 
spectrum on the cooling transition, which for a two-level 
system is of Lorentzian shape. The principle idea of EIT 
cooling is to modify this absorption spectrum such that 
the carrier transitions are eliminated. For EIT, an addi- 
tional driving field is applied as shown in Fig. [TJb) [i^] . 
If both fields have equal detuning from the respective 
qubit transition frequencies, it turns out that the qubit 
is driven into a coherent superposition of the two ground 
states \g) and |e). From this superposition, excitation 
to \a) is possible either via |e) or \g). These channels, 
however, destructively interfere, such that no absorption 
occurs even though both fields are applied. This phe- 
nomenon is known as EIT ^20i] . For suitable field param- 
eters, a cooling field absorption spectrum as shown in 
Fig. [Hd) can be achieved. Carrier excitations are sup- 
pressed since no cooling field absorption occurs at the 
corresponding frequency due to EIT. This way, ground 
state cooling is achieved even in the non-resolved regime. 

We now proceed with a quantitative analysis of the 
qubit-NAMR system. We assume that the NAMR vi- 
brates in its fundamental mode and in the plane of the 
loop. It has an effective mass Meff, length I, frequency 



u, and amplitude of the fundamental oscillation mode x. 
The quantum NAMR has Hamiltonian Hr = hvh^h, an- 
nihilation operator 6, and zero point fluctuation Xq — 
^JfijlMv. The flux qubit consists of a superconduct- 
ing loop with three Josephson junctions. Two junctions 
are identical, while the third is smaller by a factor of a. 
The two larger junctions have equal Josephson energies 
Ejx ~ Ej2 — Ej and capacitances Cji — Cj2 — Cj, 
while for the third one Ej^ — aEj and Cj3 ~ aCj, 
with a < 1. The gauge-invariant phase drops across the 
three junctions in the qubit loop are ipi, (p2 and ip^. The 
whole structure is exposed both to a constant magnetic 
field B perpendicular to the plane and to driving mi- 
crowave fields giving rise to time dependent magnetic 
fiuxs (TDMFs) $(i). Introducing (pp = {(pi + (p2)/2 
and (fim = ifi — </'2)/2 as the coordinates, the Hamil- 
tonian can be written a.s H = Hq + Hr + Hj, with Hq 
and Hr as the free Hamiltonian of the qubit and the 
NAMR, and Hj as the interaction part. The last term 
in Hj is crucial to EIT cooling because it includes an 
interaction of resonator, qubit and field. We define Cp ~ 
cos{2(pp+2TTf) and5p = sin(2(/jp-|-27r/), and obtain ifg = 
Pl/{2Mm) + Pp/{2Mp) - 2Ej cos ipm cos ifp - aEjCp as 
qubit part, and Hj = a/{l + 2a)^Blx + ^{t)]Pp + 
aEj^Blx + ^{t)]Sp + aEj^'^[Blx + ^(t)fCp/2 as the 
interaction part. Here, we have introduced the momenta 
Pm = —i^d/dipm and Pp = —ihd/dtpp, and the effective 
masses of qubit Mm = 2(7$^ and Mp = {l+2a)Mm- The 
bias / = <^e/^o where $ = 27r/$o with $0 as the flux 
quantum, and $e is the static bias flux corresponding to 
the equilibrium position a; = of the NAMR. We further 
assumed moderate TDMF such that ^ = ^{Blx + ^{t)) is 
small, and expanded corresponding trigonometric func- 
tions to second order in ^. As indicated in Fig. [I{b), 
we apply two TDMF with different frequencies tOLg, ^Le 
and amplitudes Ag, Ag. The corresponding Rabi fre- 
quencies are hV,j = aEj^Aj{a\Sp\j) /2, where we have 
dropped a small contribution from the momenta. The 
LD parameters are defined as r/^D = IVg ^ Ve\, where 
77, = BlXo^{a\Cp\j)/{a\Sp\j). 

Applying the rotating-wave, the Born-Markov and the 
LD approximation (tjlu <C 1), we obtain 

p = -^[Ho + Hi,p] 

+ ^ (7g, |5> {a\Bg) p + ^ (7e, |e) {a\B,) p 
+ if (F, \g){e\B:,) p + (r^/2, |e) (e| - \g) {g\) p 
+ [N{u) + \\^ {v/Q, h) + N{v)^ {ylQ, 6+) , (la) 
H^ = -hAg\g){g\ - hA,\e){e\ + hvb^b, (lb) 
Hi = nngBg\a){g\ + Ml^B^\a){e\ + H.c. , (Ic) 
^ (7, A)p = 7/2 {2ApA'^ - A^Ap - pAt A} , (Id) 

with Bj = I + rij{b + fe^) for j G {g,e,3}. Initially, the 
NAMR occupation is = N{v) = [exp{hiy/kBT) - 
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due to the thermal environment of temperature T. The 
detunings between the TDMF and the corresponding 
transition frequencies are Ag = ^ ^Lg and = 
LJae — ijJLe- We havc redefined the transition frequencies 
to include negligible level shifts. The decay rates are de- 
fined as in Fig. [TJb) , and the pure dephasing rate of tran- 
sition \g) o |e) is denoted by F^. In Eq. (fTa|) the second, 
third and forth terms describe the spontaneous emission 
of the flux qubit. The fifth term considers an additional 
pure dephasing. The final terms include a thermal bath, 
since NAMR operate at frequencies with non-negligible 
thermal mode excitation. Q is the NAMR quality factor. 

In the LD limit, the qubit degrees of freedom can be 
adiabatically eliminated to derive a rate equation for the 
vibrational NAMR states. For the analytical analysis, 
we neglect the decay from the excited state |e) to the 
ground state \g) and their decoherence because the rates 
F and F^ can be designed much smaller than the rates 
7g and 7e. In two-photon resonance Ag = Ag = A, 
the rate equation for the average number of phonons 
= J2'^=o ninlpln) of the vibrational number states 



(W + u/Q) (n) + 



vN{v)/Q + 5A+. 



is d/dt (n) 

Here, 5A+ = {rj^jg + r,he)p'a"' /2 + vl^pr'/2 is negh 

gible since pa''\ F and 773 are small, pa''^ and pi'"*-' are 
the steady-state population of auxiliary state \a) and |e) 
in the absence of the NAMR. We identify W = A_ - A+ 
with a net cooling rate in the zero-temperature case [tI, 
such that cooling requires W + v/Q > 0. The two tran- 
sition rates A± describing the heating and cooling exci- 
tations are given by [il = (17^ -t- fll)^/^] 



A+ = 



f]2^ 72;/2+4[l]2_j,(j,±A)]^ 



(2) 



The steady state evaluates to n 
and is minimal for Q, 



vNi/QW + A+/W 



^v{v- A) and A < 0. Then, 
the cooling rate scales as Wmax '^Vld^I^IH^'^i)- 
analyze the cooling dynamics further, we assume | A| > 
define r = fie/^^g, and distinguish the two cases of weak 
(r ^ 1) and strong cooling fields {r k, 1). We find 



nss{r > 1) 
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(4A)2 
2 
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(3a) 



(3b) 



r/2^Q|A| (4A)2 ■ 

Similar to sideband cooling, Uss has two contributions. 
The first arises from the environment and is proportional 
to iVi, while the second one is due to cooling field scat- 
tering. But there are crucial differences to other cool- 
ing approaches. Even though we operate in the non- 
resolved regime, the field scattering contribution is in- 
versely proportional A2, allowing for ground state cool- 
ing. In contrast, this is only possible in resolved sideband 
cooling 0, Sj- Also, other than in backaction cooling Q, 
no significant coherent shift of the final phonon num- 
ber occurs in FIT cooling. Finally, in FIT cooling, the 



environmental contribution to the cooling limit is sup- 
pressed by a factor of A relative to sideband cooling [7| 
and backaction cooling. For the latter comparison, we 
define an effective LD parameter rj^D = rj^yrimax, where 
V = (xq/u) d/dxujc with resonance frequency ujc of the 
cavity and rimax as the photon number in resonance Q. 

We now turn to our results based on our analytical 
analysis and on a numerical solution of the full system 
Eq. ((Ta)) 2JI . We choose the fiux qubit parameters based 
on the experimental work 25| as a = 0.7, Ej = 200GIIz, 
and Ej/ Ec=50, where Ec is the charging energy of 
the junction. We assume a bias flux / — 0.5005, and 
calculate transition frequencies oJeg ~ 27r x 4.89GIIz, 
LOag w 27r X 30.68GHz and Uae « 27r x 25.79GHz. The 
decay and pure dephasing rates are taken as the mea- 
sured values F - 27r X 2MHz and F^ - 27r x 4MHz at 
/ = 0.5005 |25j]. The linewidth of state \a) is given by 
7 = 7g + 7e. Since 7 was not measured in [25|, we as- 
sume 7 — 50F [2^ [23] . As the final phonon number is 
insensitive to 7g/7e, we choose 7g = 7e = 7/2. The LD 
parameters are rjg = 28.19BIXq^ and r/e — —OAOBIXq^ 
and r?3 = BlXo^{e\Cp\g) / {e\Sp\g) = 0.02B/Xo$, respec- 
tively. Rabi frequencies Jig ~ 7 can be achieved with 
input power < 0.5 /jW [2^. This power can be fur- 
ther reduced via a larger mutual inductance between 
qubit and driving circuit. The NAMR is taken as a dou- 
ble clamped resonator with size 25 /itmx lOOnmx 120nm, 



effective mass M, 



eff 



< 



2pg, resonance frequency 



2tt X 25MHz and quality factor Q = 5 x 10^ [Ulii,!!!]. 
We assume a magnetic field B < 3T, which is resonable 
in NbN-based qubits, to obtain rjLu = 0.0566. 

First, we study the dependence of the cooling limit 
on the temperature of the environmental bath. Cor- 
responding results are shown in Fig. for both weak 
and strong cooling fields. In the "small TV^-large rjj^uQ" 
regime, the dominant contribution to the cooling limit 
arises from the scattering of the cooling TDMF. Ne- 
glecting the environmental contribution, the steady-state 
phonon number Uss then simplifies to Uss — 72/(4A)2. 
This result is independent of r, unless r becomes large 
enough to outweigh the small ratio Ni/ {rj^j^Q). Conse- 
quently, saturates to a constant value towards smaller 
Ni. In the second "large iVi-small rj'^^jjQ" regime, the 
contribution proportional to Ni arising from the envi- 
ronment dominates the cooling limit. In Fig. [2Ia), we 
see that our approximate analytical results agree well 
with the numerical results for negligible decoherence on 
the ground state transition. The dash-dotted curves in 
Fig. [2lja) show the numerical results with decay and de- 



phasing rates as measured in [25|. With decoherence, 
field parameters slightly modified from the analytical pre- 
dictions for the decoherence-free case improve the cool- 
ing. For example, assuming TV^ = 16 corresponding to 
a temperature T — 20mK, a weak [strong] cooling field 
cools the system down to Uss = 0.65 [n^s = 0.71] in this 
case, i.e., close to the motional ground state. 
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occupation. As in the single-mode case, cooling improves 
further with decreasing F and F^. 

In summary, we studied ground state cooling of a 
nanomechanical resonator coupled to a flux qubit in the 
non-resolved regime. Efficient cooling is achieved because 
detrimental carrier excitations are suppressed by quan- 
tum interference. 



FIG. 2: (Color online) (a) Cooling limit Usa against initial 
phonon number A'^;. Dashed lines show Eq. (|3]), and solid 
lines show full numerical results with F = 0, F^ = 0. The 
dash-dotted lines show numerical results including decoher- 
ence F = O.O27, F,^ — O.O47 as measured in [25|. Results 
are shown for two different cooling field strenghts. Param- 
eters are Vl^ = \/ v(y — A) v(y — A)/2] for weak [strong] 
cooling fields and A = —87, except for the dash-dotted curve 
with r2e = which shows Ag = —2.857, = ~37, and 
Jle = 0.537. (b) Cooling limit a function of the detun- 



Eq, 



((3|, and (ii) numerical results for F = F^ = 0, both with 
= n""* := ^Jv{v-Ae)/2 and A^ = Ag. The other curves 
are numerical results for Ag = 0.99Ae, F = O.O27, and (iii) 
F^ = 0, fie = 0.58n°f', (iv) F^ = 0.047, fie = 0.64n°f* and 
(v) F^ = O.O87, fie = 0.64^°''*. 



In Fig. [21b), the dependence of the steady state phonon 
number n^s on the detuning A is shown for T — 20mK 
{Ni = 16). From our analytical results, we expect 
stronger cooling fields to be particularly effective if larger 
detunings are used, since the environmental contribution 
is suppressed by the detuning A, see Eq. ([Sb]) . There- 
fore, we focus on stronger cooling fields (r = 1). As 
before, good agreement between theory and numerical 
results is obtained for the case with negligible ground 
state decoherence. As predicted from the analytical re- 
sults, the cooling limit decreases with increasing detuning 
|A|. Note, however, that our analytical results become 
invalid for |A| > IO7 because the system then cannot 
be excited efficiently. For the numerical results including 
decoherence, we obtimize the chosen driving field param- 
eters slightly away from the parameters suggested from 
the analytical calculation. For the decoherence parame- 
ters measured in 25| shown in curve (iv), a cooling limit 
oiuss — 0.44 is achieved, corresponding to a ground state 
occupancy of 70%. As expected, higher [lower] decoher- 
ence rates lead to higher [lower] final phonon number. 
Note that in principle, F and F^ can be smaller by several 
orders of magnitude than 7 in our system [S^] , leading 
to a cooling performance as in curve (ii) of Fig. [21 

Finally, we note that our scheme can simultaneously 
cool several resonator modes For example, assuming 
F = 0.027, F0 = 0.047, Ae = -37, Ag = -2.867, and 
U,g — D.f, — 0.6947, the point of vanishing absorption 
moves between two modes vi — 7/4 and V2 = 7/2. We 
find that one can then simultaneously cool both modes 
from Ni = 16 to about Uss = 1, i.e., 50% ground state 
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